Engineered optical nonlinearity for quantum light sources 
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Many applications in optical quantum information processing benefit from careful spectral shaping 
of single-photon wave-packets. In this paper we tailor the joint spectral wave-function of photons 
created in parametric downconversion by engineering the nonlinearity profile of a poled crystal. We 
design a crystal with an approximately Gaussian nonlinearity profile and confirm successful wave- 
packet shaping by two-photon interference experiments. We numerically show how our method can 
be applied for attaining one of the currently most important goals of single-photon quantum optics, 
the creation of pure single photons without spectral correlations. 



Spontaneous parametric downconversion is a nonlinear 
optical process in which a photon from a pump laser, inci- 
dent on a nonlinear birefringent crystal, converts into two 
single photons under conservation of energy and momen- 
tum. Photon sources based on this phenomenon are an 
ubiquitous tool for quantum computation [l], quantum 
communication [2 and quantum metrology [3l [4] . They 
are also becoming increasingly important in more spe- 
cialised applications such as quantum imaging [5], quan- 
tum lithography [6] or optical coherence tomography [7] . 
As these experiments evolve, more stringent requirements 
are placed on the characteristics of the quantum state 
of emitted light. In particular, to produce high-purity 
heralded, or even near deterministic single photons, the 
spectral shape and correlations of the created photon 
pairs must be carefully engineered. 

The most common method for spectral engineering 
is filtering, however, this can lead to loss and mix- 
ing. A more sophisticated method involves shaping the 
spectrum at the source. Consider downconversion in a 
quasi-phasematched (QPM) crystal with a poling pe- 
riod A [8 , where the crystal domain is inverted when- 
ever the pump and downconversion fields acquire a phase 
mismatch Ak=2m7r/A — where m is an odd integer — 
allowing phase-matching of a wide range of wavelengths 
in different nonlinear materials. There are several meth- 
ods to control the joint spectral amplitudes of photons 
created in downconversion in a QPM crystal [H [9] . For 
example, imposing a linear chirp on the poling period A 
has been used for the generation of ultra-broad-spectrum, 
top-hat shaped photons [7] for optical coherence tomog- 
raphy. However, the currently known methods involve 
changing A, which may be incompatible with stringent 
phasematching conditions. 

In this paper, we consider type-II downconversion in a 
QPM crystal with a longitudinally non-uniform grating. 
We can synthesise photon pairs with arbitrary spectral 
amplitudes by modulating the nonlinearity profile x(^) 
of a crystal through different-order poling without chang- 
ing the phase-matching conditions. We tailor a spectral 
photon-pair amplitude with an approximately Gaussian 
profile, which is critical for high purity photon genera- 
tion, as we will discuss later. In addition, Gaussian spec- 



tra are known to be optimal for temporal mode matching 
|1QJ — a critical consideration in any experiment involving 
single photons. 

Theoretically, the two-photon component of the optical 
state is described by [11] 

^p) = dujidujsf{uji,ujs)dl{uji)al{ujs)\0)^ (1) 

where f{uji^ujs) = a{uJi-\-uJs)^{oJi^oJs) is the joint spectral 
amplitude of the created photons in the idler and signal 
modes respectively (for details, refer to Appendix [A] ). 
The spectral properties of downconverted photons can be 
manipulated via the pump envelope function a{uji + ujg) 
[12 , or as we show here, the phase matching function 
(PMF) ^{uji^ujs)' We use a monochromatic pump a{uJi-\- 
Us) = S{uJi -\- oJs — l-^p) where jip is the pump frequency, 
and directly tailor 
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which is the Fourier transform of x(^) ^^^ Ak is the 
phase mismatch. 

According to equation [2J the phase matching function 
of a standard crystal with a uniform nonlinearity pro- 
file is ^{uji^ujs) = sinc(AA:L/2). However, to generate a 
Gaussian phase matching function, we require a crystal 
with a Gaussian nonlinearity. While it is non-trivial to 
directly change the material properties, we can make use 
of higher-order poling to realise a variety of nonlinearity 
strengths. In a poled structure, the effective nonlinear- 
ity scales with the poling order m as Xeff=2x/'7r^Ti. For 
odd 771, mth order QPM can be achieved by reversing 
the direction of the poling every m coherence lengths, 
defined as Lc = A/2. Even-order QPM can be achieved 
by combining two odd orders. 

We exploited this feature to design a crystal consist- 
ing of a number of discrete sections, each with a differ- 
ent Xeff, discretely approximating the desired Gaussian 
shape, and custom-poled a 10 mm long Potassium Ti- 
tanyl Phosphate (cpKTP) crystal accordingly. Details 
on the actual design of our tailored crystal can be found 



in Appendix B Our technique decreases the overall non- 



linearity and reduces the effective length of the structure. 
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FIG. 1: a) Nonlinearity profile for tfie cpKTP crystal Xt(^) (orange line) and target Gaussian profile Xg{z) = exp( — (2;/Leff)^/7) 
(black dashed line) with effective length Lgff = 5.67mm (green dot-dashed line) and 7 ~ 0.193 (see Appendix C). b) Phase- 



matching function amplitudes and intensities (inset) for the cpKTP (orange line) compared to a ppKTP of the same effective 
length Leff (green dot-dashed line) and target Gaussian profile $g(<^z, <^s) = exp(— 7(AA:L/2)^) (black dashed line), c) Magnified 
image of part of the custom-poled KTP crystal. Vertical lines separate sections with constant effective nonlinearity, with their 
poling order ??i, length L and poling duty cycle D. d) Magnified view of the transition from poling order ?7i=l to ?7i=2 [? . 
Due to a slight mismatch between design and actual domain lengths, the crystal was shortened by a few tens of fim. on one 
side. 



Appendix C We therefore expect less photon-pair yield 



and broader bandwidths when compared to a standard 
periodically poled KTP (ppKTP) with the same length 
and phase-matching. 

Fig. flk) shows the tailored nonlinearity profile Xt{z)^ 
defined in Eq. (Bl) in Appendix [BJ together with the 
target Gaussian profile Xg(^)- The corresponding PMF, 
obtained from the inverse Fourier transform of Xt(^), is 
very similar to a Gaussian function, as shown in Fig.fit)). 
Compared to the sine-shaped phase matching function of 
a ppKTP of the same effective length as the cpKTP (5.67 
mm), the side lobes on either side of the central peak are 
significantly suppressed. This becomes even more evident 
when considering the spectral intensity (see inset). 

A microscopic image of part of the custom-poled KTP 
(cpKTP) crystal is shown in Fig. flp). One can clearly 
see the individual sections with different poling orders, 
which line up with the theoretical design almost perfectly. 
Fig. Illi) shows a magnified view of a transition between 
poling-order sections m=l and m=2. 

We tested our custom-poled crystal in a typical down- 
conversion setup, see Fig. |2k), comparing it to a 10 mm 
long ppKTP crystal (A=10.95 jam). Due to the reduced 



overall effective nonlinearity, we expect a relative photon 
pair rate of 34.4%. The measured rate (detected with- 
out a beam-splitter) was '^10 kpairs/s for the custom- 
poled and ~33 kpairs/s for the standard crystal, respec- 
tively. This corresponds to a relative yield of ~ 30.4%; 
we attribute the small reduction in efficiency to the fact 
that the custom-poled crystal, in contrast to our stan- 
dard crystal, was not anti-reflection coated. 

We verified the spectral amplitude of bi-photons cre- 
ated in the cpKTP in a two-photon interference exper- 
iment. When two indistinguishable photons mix on a 
symmetric beamsplitter, they will always be found in the 
same output port. This phenomenon was first reported 
in the landmark experiment by Hong, Ou and Mandel 
[14^ , who observed a dip in the photon-pair detection 
probability as a function of the temporal delay between 
the input photons. Theoretically, the shape of this co- 
incidence dip assumes the inverted Fourier transform of 
the absolute square of the PMF, which can be readily 
calculated from the two-photon state in Eq. (IT]), see fl4]. 

The recorded interference patterns for the cpKTP, and 
the standard ppKTP are shown in Fig. [2b), along with 
the theoretical interference patterns calculated straight 
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FIG. 2: a) Experimental scheme. The crystals (cpKTP, ppKTP) were temperature-stabilised (TEC) and pumped by a 410 nm, 
grating-stabilised diode laser. The emitted orthogonally polarised photon pairs were split at a polarising beamsplitter (PBS) 
and coupled into single-mode fibres equipped with polarisation-controllers (POL). They were then superposed at a 50/50 fiber 
beamsplitter (BS) and detected in coincidence. We obtained two-photon interference patterns by changing the delay At with a 
motorised translation stage. The only filters in use were two RG715 long-pass filters (LP), b) Two-photon interference patterns 
for the cpKTP (red circles) compared to a standard ppKTP (green diamonds). The solid lines show the theoretical values, 
calculated from the respective PMF for each crystal. The reduced chi-square values of these fits are 3.07 and 5.51, respectively. 
The dashed lines show least-square fits of a triangular pattern to the tailored crystal data and a Gaussian fitted to the normal 
crystal, with reduced chi-square values of 50.59 and 23.10, underlining the strong divergence from these shapes, c) Spatial 
quantum beating for various center- frequency detunings Acj — uji — ujs. The lines show the ideal values, calculated from the 
respective PMF. All probabilities pc for b) and c) were obtained by normalising detected pairs to twice the averaged counts 
outside the coherence length. The only free parameter for theory values was the interference visibility of ^ 95%. All error bars 
are smaller than symbol size. 



from the PMFs in Fig. lip). The bandwidth difference 
results from the different effective lengths of the two crys- 
tals. The interference pattern for the standard crystal is 
triangular, just as expected for the sine-shaped PMF [15]. 
The pattern for the custom-poled crystal departs from 
the triangular shape and indeed approximates a Gaus- 
sian. The interference visibility was ^ 95% for both crys- 
tals, confirming that the indistinguishability of the down- 
converted photons was not compromised by the crystal 
modulation. 

To further explore the underlying spectral correlations 
in the PMF, we measured spatial quantum beating pat- 
terns. We detuned the center frequencies of the downcon- 
version photons via a change in crystal temperature away 
from its optimal value for collinear, degenerate quasi- 
phasematching and again observed two-photon interfer- 
ence [15]. The results in Fig. ^) show that the custom- 
poled crystal exhibits less distinct beating, in particu- 
lar, less anti-bunching, i.e. coincidence probability values 
above the random level of 0.5. The maximum value for 
the cpKTP was 0.546 ± 0.005 compared to 0.586 ± 0.003 
for the standard crystal, a significant reduction relative 
to the base-line of 0.5. The theory values for the ppKTP 



were adopted from [15] directly. For the cpKTP we re- 
peated the calculation in [15] using the respective custom 
PMF shown Fig. [i]d). 

The observed anti-bunching occurs when the 
frequency-detuned spectral wavefunction (i.e. the 
joint spectral amplitude) of the two-photon state is 
partially anti-symmetric, which in turn reveals the 
frequency entanglement intrinsic to downconversion [15]- 
[TTl . A Gaussian spectral amplitude is always positive 
and therefore does not have anti-symmetric components, 
which explains the significantly reduced beating in the 
interference patterns of the custom-poled crystal. All 
measured interference patterns for our custom crystal 
agree exceedingly well with the patterns calculated from 
the theoretic PMF. 

One situation in which this method will be useful is 
the generation of pure heralded single photons. Even in 
a group- velocity-matched configuration with a symmetric 
joint spectral amplitude, the maximum purity of the her- 
alded single photon state is, due to the PMF sine profile 
|18j . limited to 0.81 (for ppKTP with a 788 nm pump 
[r9-21 without spectral filtering). This purity affects 
the two-photon interference visibility between heralded 
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single photons and thus ultimately the quality of opti- 
cal quantum gates or multi-photon states generated in 
post-selection. We numerically compare two designs for 
cpKTP crystals (see Appendix |e|, characterised by their 
maximum poling order {m=l and m=2)^ with a stan- 
dard ppKTP crystal following [18] . Table [l] shows the 
calculated purities for our two designs compared with 
a standard crystal. Taking 7Timin=l yields a coarse ap- 
proximation to the Gaussian function, but the purity P 
improves substantially to 0.97. A better approximation 
to the Gaussian is achieved if mmin=2, then P=0.99. For 



further detail, refer to Appendix E 



TABLE I: Numerical comparison of the purity P of heralded 
single-photons of a standard crystal and two cpKTP crystals 
of length L, in a group- velocity-matched scenario. The effec- 
tive length of both cpKTP crystals is Leff=24.2mm. 



Crystal 



L [mm] Purity P 



ppKTP 
cpKTP m 
cpKTP m 



mm 



mm 



24.2 0.81 

1 40.5 0.97 

2 41.6 0.99 



In conclusion, we demonstrated longitudinal shaping 
of single-photon wave-packets via indirect modulation of 
the nonlinearity of a crystal. Our method can be used 
to generate other spectral profiles of interest, such as a 
triangle or a top-hat, see Appendix |FJ A comb-like non- 
linearity structure, for example, would allow the direct, 
lossless generation of frequency-bin qubits [22l |23]. In 
addition, the custom-poled crystals can be used in single- 
crystal sources for polarization entanglement, such as the 
Sagnac-type confgurations reported in [24l |25] , or even 
in a single-pass scheme, combined with the related tech- 
nique of using two interlaced first-order poling periods 
for concurrent type-II downconversion of frequency non- 
degenerate photons [26]. Furthermore, we expect this 
technique to have applications in classical nonlinear op- 
tics, e.g., in second harmonic generation, similar to the 
spectral shaping techniques previously demonstrated for 
this regime [HI |9] . 

Since our technique allows manipulation of the phase- 
matching profile, it complements group- velocity match- 
ing, which is commonly achieved by controlling the ori- 
entation and width of the phasematching function. A 
Gaussian shape vastly improves the purity of heralded 
photons from SPDC. The inevitable reduction in the ef- 
fective nonlinearity is an acceptable tradeoff given that 
modern crystals have drastically reduced pump power re- 
quirements. In addition, compared to the alternative of 
spectral filtering, one can actually pump at a much higher 
power without introducing photon number mixedness. 
This allows the creation of purer multi-photon states for 
quantum information processing, e.g. Fock states with 
high photon number [18]. 

It will be worthwhile to consider nonlinearity engineer- 
ing for four- wave- mixing photon-pair sources in photonic- 



crystal fibres, where the sine-shaped phase-matching 
function has been identified as a major problem [Tf] . 
However, group-velocity matching in these materials is 
already a non-trivial task which will inevitably be fur- 
ther complicated by modulation of the non-linearity. 
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Appendix 
Appendix A: Type II SPDC Hamiltonian 

For type-II up- or down-conversion, the evolution in- 
side the crystal is governed by the operator U{t) = 

Tex.p{ — ^J^ dtH(t)) where T is the time-ordering op- 
erator and 

H{t) = A d(jOid(jOsdujpa{(jOp)e^ ^*dl{uJi)dl{(jOs) 

'L/2 
■L/2 



is the multimode Hamiltonian p!8| [28] . where Auj = 
uji-\-ujs—ojp^ L is the length of the crystal, A is a constant 
proportional to the nonlinearity and Ak{uJi,uJs,ojp) = 
kp{ujp) — ki{uji) — ks{ujs) is the phase mismatch. Typ- 
ically, the nonlinearity does not vary over the length 
of the crystal and we can rewrite the spatial integral 

/_L/2 ^^/^*^^^^~^^^*^^"^^^^ as the Fourier transform (FT) 
of a rectangular function V^ J^ Xr(^)e~^^^*^^''^^^^(i^, 

where Xr(^) = X^^H'^(^ + ^/2) — u{z — L/2)) and u 
is the Heaviside step function. Evaluating the spatial 
integral yields the result Lsmc{Ak{uJi^uJs)L/2) where 
sinc(x) = sin(x)/x, giving rise to the following form for 
the Hamiltonian 



H{t) = AL / duJiduJsdujpa{ujp)^{Ak{uJi,(jJsj(^s)) 



iAu;t-t 



xe^^^^a}(cJi)aT(cj,)+H.c. 



(A2) 



where 



-I 

^{Ak{uJi^ujsjUJp)) = sine (-Ak{uJi^(jOs ^ujp)L] (A3) 

is the phase-matching function. Note that by picking the 
spatial integration to be centered around z = 0, it is 
possible to eliminate a phase term which would normally 



be present in equation (A2). 



For a periodically poled crystal the phase mismatch, 
Akp{uJi,ujs^ojp) = Ak{uJi^uJsTOOp) — 27r/A, includes the ef- 
fect of the periodic domain inversion in QPM. The pur- 
pose of this domain inversion is to undo the rriTT — where 
m is an odd integer — phase error accumulated by the 
pump and downconverted fields, by introducing an addi- 
tional TT phase shift to Ak at the point of inversion. For 
successful QPM, the poling period therefore has to fulfil 
the condition A = m27i / Ak^ where m is the QPM order. 

Taking the first order term of the Taylor series ex- 
pansion of U{t)\{)) and evaluating the time and pump 
frequency integrals yields the two-photon state 




t 



-0) = / I duJiduJsf{i^i,uJs)a]{(^i)al(uJs)\^) , (A4) 



where /{uji^ujs) = a{uji + uJs)smc{Ak{uJi^uJs)L/2) and 
Ak{uJi,uJs) = kp{uji +cjs) - ki{uji) - ks{ujs)- 



Appendix B: Tailoring the crystal nonlinearity 

To tailor the crystal nonlinearity, we treat each crystal 
section 5 as a rectangular function with a nonlinearity 
inversely proportional to the poling order rris. The non- 
linearity profile for the custom-poled crystal is then given 

by 



21 



s-l 



Xt(^) = 2_] ^( ~ Z_l ^^r^rA — z) Xulz V^ mrTlrA 

s=l r=l r=l 

(Bl) 
where u is the Heaviside step function, rrir is the poling 
order of the rth section, n^ is the number of domains 
within the rth section and A=10.85 /im, for type-II, first- 
order QPM of 410 nm^820 nm+820 nm. The profile 
Xt(^) is plotted in Fig. flk). 

The design of the nonlinearity profile is subject to two 
constraints. First, the nonlinearity of each section is 
limited to discrete values proportional to I/tti. Larger 
values of vn provide smoother transitions between suc- 
cessive nonlinear ites, however this leads to a greatly re- 
duced photon creation rate. Second, the width of each 
section must be an integer number of mK/2 and a mini- 
mum of 2?7iA, therefore, larger values of m may demand 
prohibitively long sections. The ratio between positively 
and negatively poled regions — known as the duty cycle 
D = l/mA where / is the length over which the sign of 
the nonlinear coefficient remains constant — was chosen 
to be 50% for odd values of m and as close as possible to 
50% for even values, as is shown in Fig. [2fc). 

While the model presented here is not strictly valid due 
to the small number of domains within each section of 
Xt(^), modeling at the domain level. Appendix [DJ shows 
good agreement with the basic model in the frequency 
range over which the detectors are sensitive. The edges 
of the spectral response function of the detectors used in 
this experiment lie in the region of negligible amplitude 
and, as opposed to spectral filters, will not give rise to 
the photon- number mixing described in \TS\ . 



Appendix C: Gaussian approximation to sine 

function 



To determine the exact shape of the target Gaussian 
function for the nonlinearity profile, we match the width 
of the desired Gaussian PMF with the sine PMF of the 
form, sinc(A/cLeff/2), that would be generated by a stan- 
dard crystal. The appropriate function is ^Q{uJi^uJs) = 
exp(— 7(A/cLeff/2)^) where the parameter 7 ^ 0.193 is 
derived from matching the FWHM of the two functions. 
We refer to Leff as the effective length, as it does not cor- 
respond to the actual length of the final Gaussian shaped 
crystal, but rather the length of the hypothetical stan- 
dard crystal. 



Appendix D: Detailed Model 

We modeled each section of the crystal as having a 
nonlinearity inversely proportional to the poling order 
m. This approximation is only valid for a large number 
of domains in each section. Here we calculate the PMF by 
explicitly considering the contribution from each domain. 
The nonlinearity profile Xd(^) will consist of domains of 
nonlinear coefficients ±x^^\ with sign changes occurring 
at positions corresponding to the poling order and duty 
cycle (this is the case for the entire length of a typical 
periodically poled crystal). For example, in the section 
corresponding to ?7i=3, where the duty cycle is 50%, the 
sign changes every 3A, while for 772=6, where the duty 
cycle is ~ 41.6%, the sign changes from "+" to "— " after 
5 A and back again after 7A. The resulting PMF takes 
the form 



$„(Afcp) = x('^E^^(«"' 



— iAkpZj —iAkpZj-i 



), (Dl) 



where Sj is the sign of the jth domain and Zj-i — Zj is 
the width of each domain. Figure [3] shows that as A/cp 
departs from 0, the two models begin to deviate. 
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FIG. 3: Phase matching functions generated from the basic 
model (black solid line) and the detailed model (light red line). 
The inset shows a magnified portion of the PMFs, detailing 
the deviation between the models. 
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FIG. 4: a) Nonlinearity profiles and b) corresponding PMFs for: m — 1 tailored crystal (solid black line); m — 2 tailored 
crystal (dot-dashed black line); and an ideal crystal with a Gaussian profile (dashed red line). The corresponding sine PMF 
(thin dotted line) has been included for comparison. 



However, as we discuss in the main text, there is very 
good agreement between the basic and detailed models 
in the region of interest, i.e. around A/cp = 0. 



Appendix E: Crystal design for separable joint 

spectral amplitude 

Downconverted single-photons have strong spectral 
correlations which result in the degraded purity of a 
heralded state. A growing effort in engineering pulsed 
downconversion sources to produce spectrally decorre- 
lated photons includes manipulating the crystal length, 
material, bandwidth and central frequency [20", '29'-'36^ as 
well as filtering the pump field, prior to down-conversion, 
using an optical cavity [Si]. Another method imposes 
group- velocity matching to limit these unwanted correla- 
tions [191 [21]. However, the presence of side lobes — which 
arise from the sine shape of the PMF — still calls for some 
level of filtering in order to achieve high-purity single 
photons. Spectral filtering is undesirable because it low- 
ers the overall single-photon production rate as well as 
introducing photon-number mixedness which limits the 
allowable pump intensity [18]. 

We show that, in combination with our method 
of modulating the crystal nonlinearity, group velocity 
matching can be used to create high-purity single photon 
states without the use of spectral filtering. Setting the 
relationship between the group velocities of the three in- 
teracting fields such that k'^ = {kg-\-k'^)/2 and picking the 

length of the crystal to be L = \/S/'ycrp{ks — k[Y (where 

cFp is the pump width in s~^ and 7 is defined below) gener- 
ates a joint spectral amplitude (JSA), where both signal 
and idler modes have equal bandwidths. For a type-II pp- 
KTP crystal, this corresponds to a crystal length I/=24.2 
mm and a periodicity of A=68.4 /im, pumped with a 788 
nm laser with a 0.7nm FWHM which down converts to 
1576 nm in the signal and idler modes. However, for a 
standard crystal of constant nonlinearity, this will not 
result in completely pure states being generated, due to 
the side lobes in the sine function. 



To eliminate the side lobes, we want to generate a 
Gaussian PMF, ^G(^i,^s)=exp(— 7(A/cL/2)^), whose 
FWHM matches that of the PMF generated by a 
standard crystal, <l>(co'i,co's)=sinc(A/cL/2). Substituting 
^Q{(jJi^uJs) as the PMF, we can now write the JSA as 



f{uji,ujs) oc exp (- 



{ujj +cjg - 2/i)' 



2 t2 



X exp 



7Afc^L 
4 



which, due to the rotational symmetry of the two- 
dimensional Gaussian function, is separable, i.e. 

f{uji,ujs) = g{uJi)h{uJs). 

We numerically compare two designs for cpKTP crys- 
tals. Fig. [4] with a standard ppKTP crystal following 
[TS] , The results are summarised in Table I in the main 

text. 
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FIG. 5: a) Examples of tailored nonlinearity profiles, and 
corresponding phasematching functions, for a) a triangular 
phase-matching function and b) a top-hat function. The 
dashed red line shows the target functions and the black line 
the results of the discrete approximation. 
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Appendix F: Non-gaussian tAvo-photon spectra 

The method introduced in this paper can be apphed 
to the generation of almost arbitrarily shaped PMFs. As 
described above, the nonlinearity profile of the crystal 
should be tailored to the Fourier transform of the desired 



PMF. Figure [5] show examples of triangular and square 
shaped phase matching functions, as well as the required 
nonlinearity profile. Negative values of the nonlinearity, 
which are required to generate the square shape, can be 
implemented by inverting the relevant domain. 
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